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Abstract 

Background: In this paper we propose a chemical physics mechanism for the 
initiation of the glycolytic switch commonly known as the Warburg hypothesis, 
whereby glycolytic activity terminating in lactate continues even in well-oxygenated 
cells. We show that this may result in cancer via mitotic failure, recasting the current 
conception of the Warburg effect as a metabolic dysregulation consequent to 
cancer, to a biophysical defect that may contribute to cancer initiation. 

Model: Our model is based on analogs of thermodynamic concepts that tie non- 
equilibrium fluid dynamics ultimately to metabolic imbalance, disrupted microtubule 
dynamics, and finally, genomic instability, from which cancers can arise. Specifically, 
we discuss how an analog of non-equilibrium Rayleigh-Benard convection can result 
in glycolytic oscillations and cause a cell to become locked into a higher-entropy 
state characteristic of cancer. 

Conclusions: A quantitative model is presented that attributes the well-known 
Warburg effect to a biophysical mechanism driven by a convective disturbance in 
the cell. Contrary to current understanding, this effect may precipitate cancer 
development, rather than follow from it, providing new insights into carcinogenesis, 
cancer treatment, and prevention. 

Keywords: Warburg effect, Rayleigh-Benard convection, Aneuploidy, Glycolysis, 
Genomic instability 



Background 

The metabolic shift from aerobic to anaerobic glucose biochemical energy processing 
by cells is strongly correlated with the transition to cancer, or as some have come to 
characterize the process, a reversion to a more primitive and competitive level of cellu- 
lar existence (Warburg [1], Szent-Gyorgyi [2]), which may still possess some rudimen- 
tary cooperative elements e.g. early metazoans (Davies and Lineweaver [3]). Our focus 
in this manuscript is to develop a molecular physics model based on non-equilibrium 
thermodynamics to quantitatively describe that process. By better understanding this 
transition we should be able to not only address cancer more effectively but also other 
metabolic diseases including mitochondrial diseases (e.g. [4]) and diseases of proton 
pumps (e.g. [5]). This modeling approach may also shed some light on the relationship 
between the Warburg effect for cancer and the so-called inverse-Warburg effect [6] for 
neurological diseases e.g. Alzheimer's disease [7]. 

O© 2013 Rietman et al.; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative 
BiolVlGCl C6ntTcll Commons Attribution License (http://creativecommons.Org/licenses/by/2.0), which permits unrestricted use, distribution, and 
reproduction in any medium, provided the original work is properly cited. 
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The living cell is an extremely complex molecular network of tens of thousands of 
different types of molecules from ionic species and small molecules to large polymers 
and polymer networks. Naturally, the number of large, multi-nanometer size polymers 
is not as high as the smaller molecules, but a shift in their numbers can easily result in 
the emergence of a disease state. For example, variation in the number of genes 
expressed due to changes in the genome (aneuploidy), can indicate cancer. Many of 
these large polymers are enzymes, or chemical reaction catalysts. Typical molecular re- 
actions in the cell are represented by the relation: 



E + S^ES^>E + P 

kr 



(i) 



where the symbols E, S, ES, P represent the enzyme, substrate, enzyme-substrate com- 
plex, and reaction product, respectively. The coefficients labeled by ks represent for- 
ward, reverse and enzyme-substrate decomposition rate constants as indicated by their 
subscripts. Obviously, if there is a huge abundance of S and limited amount of E, the 
reaction is rate limited by the concentration of E. But if the cell is malfunctioning and 
producing an excess of E when an abundance of S is present, then the cell will increase 
the P concentration by massively parallel reactions. This is governed by the well-known 
Michaelis-Menten reaction kinetics of saturable chemical reactions [8]. This effect is 
also described by the Le Chatelier principle of reaction dynamics [9]. It states that 
chemical reactions move forward or backward so as to reduce excesses in the quantity 
of reactants or products, respectively, introduced into in the reaction vicinity. 

Given the reaction dynamics afforded by Le Chatelier s principle, it is possible to im- 
agine that an external concentration of, for example, glucose allowed entry into a cell 
would ripple through the molecular network of a cell and produce an excess of the ap- 
propriate glycolytic enzymes and other molecules associated with anaerobic processing. 
This would in turn produce an excess of lactic acid and a shift in the hydrogen ion con- 
centration. Neither of these causes will necessarily induce a cell to transition to a 
higher entropy state (disorganization) of a cancer cell [10], but we argue that these re- 
action processes can disrupt the mitochondria and/or the cytoskeleton in part via the 
microtubule growth rate which is dependent on the pH value [11], a parameter known 
to change between normal and cancer cells. In particular, in cancer cells, intracellular 
pH is alkalized [12-15], and extracellular pH is acidified [12,14,15]. These biophysical 
changes in turn could contribute to mitotic failure and in rare cases aneuploidy - pro- 
ducing a stable cell exhibiting a glycolytic shift, i.e. the Warburg effect. Our hypothesis 
is schematically outlined in Figure 1. 



The model 

We argue that non-equilibrium thermodynamics is a driving force for the ultimate 
transition of a cell from the normal to the cancer state, and that the first step in the 
process is an increased operation of the cytoplasmic glucose processing machinery due 
to a non-equilibrium mechanism analogous to Rayleigh-Benard convection (Figure 1) 
[16,17]. Our proposed scheme that describes the development of the cancer phenotype 
at a cellular level is essentially a biophysics model, which should clearly be considered 
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Healthy Cell to Aneuploidy Transition: 
Initiation of Metazoan 1.0 and Stable Warburg Effect 




y 




Figure 1 Diagram of proposed transition from healthy to aneuploidy cell. An excess of external 
glucose causes a series of phase transitions starting with glycolytic oscillators which induce a pH imbalance 
precipitating in mitochondrial membrane breakage followed by microtubule depolymerization and 
mitotic failure. 



in parallel with more biochemical schemes [18-22]. We begin by reviewing some rele- 
vant physics of nonlinear dynamics and thermodynamics. 

Rayleigh-Benard Convection 

To appreciate the relevant dynamics, we can consider a thin layer of oil, of thickness d, 
lying on top of a uniformly heated metal plate. On top of the oil is another plate of 
high heat conductivity. The apparatus has non-thermally conducting bounding walls to 
keep the oil confined. We thus have a thermal gradient from a hot zone to a cooler 
zone. Initially the heat transfer is strictly via conduction from the hot surface to the 
cooler surface. The upper plate is kept at temperature T 0 and the lower plate is kept at 
T 0 + AT (AT > 0). As long as the temperature difference remains small, frictional forces 
in the fluid, due to viscosity, will keep convective motion near zero. Nonetheless, the 
conduction results in a small upward vertical vector from the bottom layer, and due to 
gravity a small downward vertical vector from the top layer. As a result of this, there 
will be a phase transition leading to a thermal convection mode of transport and an or- 
dered state. The characteristic time for the displacement of the small regions of the 
fluid is given by: 

TRB " pgad{AT) (2) 

where ^is the dynamic viscosity of the fluid, a is the thermal expansion coefficient of 
the fluid, p is the mean fluid density, and g is the gravitational acceleration. The condi- 
tion for sustained convection is given by: 

tefm*.. ( 3 ) 

where D T is the thermal diffusion coefficient and the dimensionless constant R RB is 
known as the Rayleigh number. In this case we have used the subscript RB to remind 
us that these are the Rayleigh-Benard relations. The Rayleigh number is related to the 
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temperature differential Above this convection threshold, the fluid undergoes a phase 
transition to an ordered phase - an example of a so-called symmetry breaking 
phenomenon characteristic of many critical phenomena in physics and other fields 
[23]. What is special about this example is that it is a non-equilibrium effect driven by 
the thermal gradient across the sample. At this point one can observe convection pat- 
terns known as Rayleigh-Benard convection rolls shown in Figure 2. These patterns are 
an observable manifestation of this non-equilibrium phase transition with an associated 
symmetry-breaking effect being the spatial periodicity of the rolls and a characteristic 
time scale representing an oscillation period. 

Davies, Demetrius and Tuszynski [10] provided a qualitative argument that cancer is 
a dynamical phase transition. Here, we argue, by analogy, that a similar phenomenon 
to Rayleigh-Benard (R-B) convection, a symmetry breaking non-equilibrium phase tran- 
sition could occur for cells grown with an excess glucose gradient between the outside 
of the cell and the inside of the cell, or with a critical change in the internal glycolytic 
enzyme distribution. This could represent a pre-cancerous transformation in terms of 
biological consequences of this physical process. Analogous to the R-B characteristic 
time relation above, we may now propose that for a cell: 



0a(AC) 



(4) 



where AC is the change in concentration, and using v to represent the kinematic vis- 
cosity which has the units cm 2 sec" 1 . Water has a viscosity of about 1 cP (cgs units) at 
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Figure 2 Example of Rayleigh-Benard convection rolls. (Based on GNU Matlab code of Parmigiani et al. 
[24]). 
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ambient conditions. We will assume the viscosity inside a cell is about 1000 times 
higher (e.g. 10 P, Kalwarczyk et al. [25]). Instead of the force of gravity that has no ap- 
preciable effect on living cells at this scale, we will use energy density with units erg g" 1 
(with dimension cm 2 /sec 2 , instead of gravity's cm/sec 2 ), and represent it by the symbol, 
O. A similar parameter, energy-density-rate (erg g" 1 sec" 1 ) has been described 
by Chaisson [26] who discusses at length the application of energy-density-rate 
for multiple phenomena in the universe including evolution of life and evolution 
of stars. The parameter O is a universal parameter used to describe energy and 
the rate of energy flow; in other words, a specific metabolic rate, which repre- 
sents the rate of energy transduction per unit mass. It is essentially a measure 
of the complexity of an open system and its rate of free energy utilization. In 
the case of small resting animals the value is on the order of 10,000 erg sec" 1 
g" 1 . For a single cell it is on the order of a few hundred erg sec" 1 g" 1 . Obviously, it is a 
nonlinear function due to its relationship to metabolism. As has been well-known for 
many decades, metabolic rates scale with a fractional exponent in the range of 0.5 to 1.0, 
typically 0.75 [27]. For specific metabolic rates (metabolic rates per mass unit), these expo- 
nents become in the range from -0.5 to 0.0 [28]. Finally, we use a, which is the analog of 
the thermal expansion coefficient and therefore in this case the inverse of concentration, 
and we judiciously set its value to be 10" 4 . Keeping in mind it has the same units as inverse 
concentration but when multiplied with the concentration, as shown in Equation (4) the 
units cancel. We thus have the following units for Equation (4): 



T 



L_ ML 2 
T ~r 2 ~ 



1 _L 2 T 2 M 
M ML 2 



(5) 



Most of the remainder of the paper will be devoted to present arguments and ana- 
lyses aimed at supporting the use of Equation (4) in the context of the Warburg 
effect. 

As can be seen in Figure 3, an increase in the external concentration results 
in the time constant falling more rapidly and larger values of the energy-density 
also result in smaller time constants. As the concentration increases, the value for O 
decreases. As the cell transitions from oxidative phosphorylation to glycolytic processing, 
the complexity of the molecular network transitions from mitochondrial to substrate level 




T(S) 



external concentration C (mM) external concentration C (mM) 

Figure 3 Data from Table "Range of Glycolytic Oscillations in Yeast Extract" in Hess et al. [29] with 
~10% error bars. Curves showing oscillation time as a function of external concentration using Eq [4] and 
setting a = 10" 4 . 
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processing. The first primitive cells were anaerobic and had a lower O than later cells with 
mitochondria. Therefore, the model suggests that mitochondria are not being heavily 
exploited at this stage in the glycolytic oscillation experiment described by Hess et al. [29] . 

Glycolytic oscillators 

Several in vitro and in vivo experiments have already demonstrated that an increase in 
extracellular glucose [30,31] or lactic acid [32] concentration can lead to increased 
tumorigenesis. Further, hypoxia has been linked to causing the metabolic shift to gly- 
colysis, and to cause malignant progression, with much experimental evidence 
[33-38]. 

Goldbeter [39] discusses at length glycolytic oscillations, and uses a series of 
differential equations to model the dynamics for some of the chemical species 
involved in the reactions. The observed oscillations come about as a result of 
the fact that there is a delay in manufacture of the intermediate nicotinamide 
adenine dinucleotide (NADH) and the fact that there is a finite number of mo- 
lecular components available for the glycolytic processing. The cycle time has 
been observed in S. carsbegensis to be about 5 minutes using fluorescence of 
the glycolytic intermediate, NADH [40]. More interestingly, Hess et al. [29] 
measured the oscillation frequency for different doses of fructose or glucose as 
input. As shown in Figure 3, as the concentration of the glucose (or fructose) 
increases, the period of oscillation decreases. This can only come about from in- 
creased concentration in the number of enzymes to participate in the reactions. 
As the number of enzymes and the number of molecular components increase, 
so does the entropy because the increase in numbers enables more ways to dis- 
sipate free energy, in this case represented as chemical potential of a glucose 
gradient [41]. As the concentration increases the cycle time decreases indicating 
a more efficient processing of glucose per time unit. 

More recently, Aromolaran et al. [42] describe, and show experimentally, 
glycolytically generated adenosine triphosphate (ATP) and Ca 2+ waves propagat- 
ing through a cell from application of glycolytic inhibitors focally injected from 
a glass pipette with a 1.5 micron diameter tip. The authors discover that 
glycolytically generated ATP is likely a key modulator of Ca 2+ homeostasis. Of 
course, this has direct effect on the permeability of the mitochondrial wall, as 
shown by Yang et al. [43] who describe glycolytic oscillation depolarizing the 
mitochondrial membrane. The authors describe a model based on the logistic 
function showing there is a region where the oscillations are "too rapid for ob- 
servation." Though they do not use the term chaos', this is likely a chaotic state 
observed in the logistic functions and other "chaotic dynamical functions." 

In the case of Rayleigh-Benard convection rolls, as shown in Figure 2, the rolls can 
be modeled with a sine-circle map, for example, 0 i+1 =f(6i) where the function is peri- 
odic in the angle [44,45]. There are some theoretical arguments for the glucose oscilla- 
tors being embedded in the cell membrane. Demetrius et al. [27] argue that the 
enzymes' concentration would oscillate due to periodicities in the redox potential and 
the result could be modeled as harmonic oscillators. Further, Tyner et al. [46] measure 
electrical gradients in the cell, and we show that a relevant protein, glyceraldehyde 3- 
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phosphate dehydrogenase (GAPDH) associated with glucose processing in the cyto- 
plasm accumulate at the membrane thus showing experimental support for our hy- 
pothesis that molecular oscillators accumulate at boundaries. Lastly, Pokorny [47] 
suggests Duffing oscillators [48] as a potential model of oscillatory states of a cell. In 
the following subsections we first provide some experimental validation for the differ- 
ences in GAPDH localization. Then we follow that with some simplistic analytical 
modeling to show phase locking of oscillators potentially resulting chaos and disrup- 
tions of mitochondria. 



Experimental validation for GAPDH distribution 

These predictions of GAPDH distribution were tested by using antibodies against 
GAPDH in normal human mammary epithelial cells (HMECs) under normal culture 
conditions, purchased from Life technologies. 

Briefly, HMEC cells were grown on coverslips and incubated at 37 C, 5% C0 2 
until they reached 50% confluency. They were then fixed with 4% paraformalde- 
hyde, permeabilized with 0.1% Triton and blocked for 1 hour in 1%BSA/1X PBS. 
The cells were then incubated with 1° Antibody (GAPDH from Santa Cruz #sc- 
25778) at 1:100 for 3 hours at room temperature. After 3 washes with IX PBS, 
cells were then incubated with 2° Antibody (Alexa Fluor 647 anti-rabbit from 
Invitrogen A21244) for 2 hours at room temperature. After 3 washes in IX PBS, 
some cells were incubated for 5 minutes with a dilute solution (300 nm in PBS) of 
with 4',6-diamidino-2-phenylindole (DAPI) which is a fluorescent stain that binds 
strongly to A-T rich regions in DNA and used to demarcate the nucleus. After 3 
washes with IX PBS, coverslips were mounted on slides with Prolong Gold 
(Invitrogen P36934). Cells were then imaged using the Leica SP5 Confocal 
Microscope. As shown in Figure 4, intracytoplasmic GAPDH was present in 
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Figure 4 Intracellular localization of GAPDH in HMEC cells on glass slides with and without glucose. 
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highest concentrations in the cell periphery and virtually absent from the peri- 
nuclear regions of the cytoplasm. 

In some experiments HMEC cells were seeded on glass slides or Cytoo Chips at a 
concentration of 30 K cells/ml and incubated overnight at 37 C The next day, cells 
were glucose starved for 1 hour. Control cells were incubated in complete media. After 
1 hour, cells were washed in IX PBS and fixed for 10 minutes in 4% Paraformaldehyde, 
permeabilized in 0.1% Triton for 10 min, and blocked in 1%BS A/PBS for 1 hour. For 
the immunostaining procedure, cells were incubated in a 1:100 dilution of GAPDH at 
room temperature for 1 hour, followed by incubation with the secondary rabbit anti- 
body, Alexa Fluor 594, at 1:500 for 30 minutes at room temperature. The cells were im- 
aged with a Zeiss Confocal Microscope. Figure 4 shows intracellular localization of 
GAPDH in HMEC cells on glass slides with and without glucose while Figure 5 shows 
intracellular localization of GAPDH in HMEC cells on Cytoo Chips with and without 
glucose. 

These clear experimental results support our conjecture that there is an increase in 
the glycolytic processing as a result of an external glucose gradient. We now describe 
some analytical models outlining the phase locking of the glucose oscillators and the 
implications. 



Sine circle Map modeling 

As suggested in the above experimental work the glucose oscillators likely accumulate at 
specific regions in the cell to make maximum use of the large glucose gradient. To model 
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this phenomena, we will use an array of sine-circle maps in a grid, called a coupled map 
lattice, with a global coupling parameter representing the external glucose concentration. 
This global coupling will result in increases in oscillation as the coupling constant in- 
creases. The sine-circle map, represented as a difference equation, is given by 



ft 



-6 t + — sin(27T(9,) 
In 



(6) 



As k increases from 0.5 to 4.5 the output of the function transitions from a fixed 
point to 2-cycle, 4-cycle and finally chaos. This is a typical bifurcation diagram similar 
to that produced by the logistic function [49]. If we replace the 2 in the argument of 
the sine function with 3 or 4 it simply moves the bifurcation points on the 6-k curve. 
This is shown in Figure 6. 

There has been research showing a quasi-periodic route to chaos in biochemical 
systems [21,39,50,51]. The biochemical system, in these cases, glycolysis, starts at 
some given frequency and as external glucose is pulsed, starting at a frequency 
similar to the glycolysis value, the biochemical reactions begin to phase lock with 
the external pulses and first results in increased amplitude followed by period 
doubling as the external glucose pulse frequency increases. Finally, at a stochastic 
external glucose pulse frequency, the observed internal pulses are only about 3-4 
X higher than the original glucose oscillations, but the actual sine wave appears as 
if there is a frequency and/or amplitude modulation. The oscillations observed in 
the cell extracts are not chaotic. Analytical models of these biochemical systems 
[52,53], based on differential equations, support a periodic-doubling route to chaos, 
though in the real biochemical systems chaos is rarely observed. This is likely due 
to the fact that a cell is not just "a bag of chemicals" but really a hybrid between a 
chemical network and molecular nanomachines. The equations of dynamics 
also apply to mechanical systems, yet real machines are not capable of being 
driven into chaotic state, because of friction and physical constraints in the 



0.9 
0.8 
0.7 
0.6 

n 

J 0.5 
0.4 
0.3 
0.2 
0.1 



n; r Hji 

hi ; ill h 



III l I | ill Z. T II M 

: i jl j 
m li i! 



i. 



i! 




0.5 



1.5 



3.5 



2 2.5 3 

K 
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machine. Similarly, nanomachines are not likely able to be driven into chaotic 
states because of secondary bonding effects, e.g. van der Waals forces, and steric 
hindrance and potential energy surfaces [54]. 

In the above, cell-free extract, experiments, the investigators monitored NADH 
to observe the oscillations in the range of 0.002 to 0.005 Hz. Yang, et al. [43] 
did experiments on tissue models of rabbit ventricular myocytes and observed 
adenosine diphosphate (ADP) oscillations. Of course, when combining these ex- 
perimental observations it must be remembered that there are 2 NADH mole- 
cules and 4 molecules of ADP per glucose molecule in the glycolysis reaction. 
Yang et al. [43] observed oscillation increases from about 0.02 to 0.067 Hz; but 
like the cell-free experiments, they did not observe chaos. Since chaos has been 
observed in real-world chemical systems [55], the lack of observed chaos in 
these biochemical and tissue experiments can be explained as follows. 

In the case of tissue-based experiments, the oscillators are operating at a fixed fre- 
quency all driven by chemical kinetics and Le Chateliers principle. The observation of 
increase in frequency is likely the mean-field effect from an increase in the numbers of 
oscillators and this would also account for the observed effect on decrease of energy- 
density-rate as the concentration increases. This means the cell is adapting to the ex- 
cess external glucose by producing more glycolysis oscillators. As will be shown shortly, 
these oscillators can phase lock with each other and produce oscillations at a frequency 
about 2 or 3 times higher. 

In the case of the cell-free extract, it is not likely that more oscillator- 
components are being produced on demand, so Le Chatelier s principle will not be 
modulating the overall molecular network. Instead, the existing molecular compo- 
nents for oscillator construction are fixed, and more in situ oscillators may form 
because of the excess glucose. Again these oscillators can phase lock and produce 
the observed frequencies. 

We can describe the phase locking with a coupled map lattice. Since the glu- 
cose oscillators are modeled here as a sine-circle map, we build our coupled 
map lattice from these. The definitive reference on coupled map lattices is by 
Kaneko [56]. Coupled map lattices (hereinafter CML and not to be confused 
with chronic myelogenous leukemia), are lattice models with, usually, difference 
equation mapping relations in the cells comprising the lattice. The cell updates 
are given by x(i) <—f[x(i)). And to include "diffusion" or coupling between the 
cells one typically modifies the update equation as 

(i)^(l-e)f(x n (i))+-[f(x n (i+l) +f(x n (i-l))] (7) 

This is a one-dimensional CML, where the left and right neighbor of cell i 
are coupled to cell i. We use the sine circle map as the main function / 

d t+l ^e t +^sm{2ne t ) (8) 
In 

Rather than use a global coupling parameter, e, as is usually done, we assume 
a self-regulatory threshold dynamics. The adaptive mechanism is triggered when 
a glycolytic oscillator exceeds a critical threshold x% and excess is passed on to 
its neighbor. As observed in spin-glasses [57], we assume a symmetry breaking 
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effect, so that only one neighbor actually receives the excess and which neigh- 
bor, (e.g. left or right) is preserved throughout the dynamical update. Our algo- 
rithm for a one-dimensional array is thus: 

1. for each oscillator compute the current state 

/ c 

x(i)<—x(i) — sin(27T#(/)) 
2tt 

2. check threshold 

if x(i) > 

deltaif) = x(i)-x* 

else 

delta(i) = 0 

3. send excess to neighbor 

x(i + l)<—x(i + 1) + delta(i) 

4. reset delta array 

delta(i) = 0 

Hence, the adaption is triggered when x(i) > x*. This causes unidirectional 
transport S(i) = x(i) - x* ; x(i + 1) = x(i + 1) + 8(i). This algorithm has been shown 
to be capable of universal computation [58]. It does, however require careful 
tuning of the threshold and bifurcation parameters. For example, since the 
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values of the sine function can reach 1.0 and if x* = 0.15 and k = 1.0, then the 
map can blow up. The full phase diagram for x* and k is given in Figure 7. As 
expected for any chaotic attractor, there are regions of fixed point, complex os- 
cillations (e.g. 2-cycle, 4-cycle) and regions we label as undefined because one 
or more oscillators blew up. 

Kaneko [56] discusses many options of coupling the oscillators. In all cases the phase 
diagram contains regions of fixed points, oscillations, and undefined or chaos. This is 
to be expected from the bifurcation of the function as shown in Figure 6. The hard 
boundary in our phase diagram at k = 2 is to be expected as shown in Figure 7. The fact 
that the oscillations begin prior to 2.5 is due to the fact that the thresholding and signal 
transfer changes the dynamics. The regions in the phase diagram labeled n are for n- 
cycle, or strange attractor [44]. 

In mapping this phase diagram to glycolytic oscillations we would not expect all 
threshold values to be valid. If we set the threshold to 0.5 then the bifurcation param- 
eter represents the glucose dosage and we have 2 cycles over a small range until the 4- 
cycles followed by continued increase in the glucose (bifurcation parameter) results in 
n-cycles. A 4-cycle in this phase space would look like a doubled 2-cycle, amplitude 
modulated glycolytic oscillator; and an n-cycle would look like an amplitude or fre- 
quency modulated 2-cycle or 4-cycle. These types of modulations have been observed 
by Hess et al. [29] and von Klitzing and Betz [59]. 

Mitochondrial transition 

We now turn our attention to the mitochondria and explore the possible implications 
from glycolytic oscillations on mitochondrial stability. It is known that the glycolytic os- 
cillations result in similar pH oscillations (e.g. Hess et al. [29]). But these oscillations 
are about tt/2 out of phase. This has implications on the polarization of the mitochon- 
drial membrane. The cell contains buffer mechanisms to minimize pH imbalance, but 
too much polarization of the mitochondrial membrane will cause the mitochondria to 
break down. The phase lag in pH can lead to a potential problem. 

Galante et al. [20] showed that an increase of protons (decrease in pH) to an assay of 
heart mitochondria results in a decrease in the Michaelis-Menten rate constant K mt 
and an increase in the forward velocity, Vf, of the reaction 

NADH + NADPoNAD + NADPH (9) 

The ratio of VflK m shows an actual phase transition at pH 7.5. The authors did not 
discuss the significance of this phase transition. Later work by Aromolaran et al. [42] 
showed waves of Ca 2+ ions traversing the cell as a result of a localized ATP perturb- 
ation. These waves are able to traverse the entire cell within 30 seconds - far faster 
than diffusion. Other work by Ramanujand and Herman [60] show a nonlinear scaling 
of glucose metabolism in normal and cancer cells, where the scaling exponent is dif- 
ferent for both types of cells. This is analogous to our observed O variation as a func- 
tion of glucose. Lastly, Aon et al. [61] describe experiments on percolation and 
criticality [62] in mitochondrial networks of a cell. The authors used a local perturb- 
ation induced by a two-photon laser excitation. They observed a cell-wide transition to 
take place within 4 seconds resulting in depolarization of the majority of the mitochondria 
in the cell. 
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This depolarization of the mitochondria membrane also accompanies an ATP deple- 
tion. This can in turn effect the following reactions [63] : 

• D-glucose + ATP — > D-glucose-6-phosphate + ADP 

• D-fructose-6-phosphate + ATP — > D-fructose-l,6-bisphosphate + ADP 

• 1,3-bisphosphoglycerate + ADP — > 3-phosphogly cerate + ATP 

• phosphoenolpyruvate + ADP — > pyruvate + ATP 

with the following heats of formation, respectively: 

AG = -16.7 kj/mole, AG = -14.2 kj/mole, 
AG = -18.8 kj /mole, AG = -31 .4 k) /mole 

ATP is needed to maintain Ca 2+ . A lag in production of ATP, as the above reactions 
compete with other reactants and products in the overall molecular network, could in- 
duce changes in the cytoskeleton via pH effects on the growth dynamics of the 
microtubules. 

The dynamics of the entire system of substrate metabolism of glucose is beyond the 
scope of this paper and has already been reported in the literature [64-67]. But even a 
small system of two coupled components, where one is phase lagged with respect to 
the other can give rise to very complex dynamics. For example, the system given by: 

dx/dt = ax — bxy , . 

dy/dt = — cy + dxy 

can result in essentially positive feedback. We show this in Figure 8. Gehrmann et al. 
[68] model the complete glycolysis reaction system and shown that the enzyme, phos- 
phofructokinase (PFK) is a key component for not only stabilizing ATP but key because 
it can also result in positive feedback. This in turn can result in a runaway oscillation 
system similar to that shown in Figure 6. The authors also show chaos, but that is likely 
an artifact of the modeling since real chemo-mechanical systems comprising nanoscale 
machines operating on molecules could not possibly result in chaos. More likely the 
system would break. In this case apoptosis would be the result, not chaos. 
It is interesting to note that our analysis suggests that ATP oscillations could effect 

2phosphoenolpyruvate + ADP— > pyruvate + ATP 

amongst other reactions. Using Le Chatelier s principle we can argue that forcing this 
reaction in the reverse direction by using 3-bromopyruvate will deplete the ATP and 
thus induce the cancer cell to enter into an apoptotic state. This type of behavior has 
in fact been reported by Ko et al. [69] and Mathupala et al. [70], who reported signifi- 
cant reduction in tumor volume in mice by treatment with 3-bromopyruvate. 

Cytoskeleton disruption 

The metabolic shift to glycolysis leads to acidosis, which subsequently results in an 
acidic extracellular pH [12]. Gillies has documented an acidic extracellular environment 
in numerous tumors [12]. Several studies link acidosis to genomic instability. Morita 
et al. show low pH leads to sister-chromatin exchanges and chromosomal aberrations 
[71,72]. Brusick, Cifone, and Cipollaro report that low pH environments (~pH 6.5) 
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caused genomic toxicity [73-75]. Several studies link the tumor environment, with hyp- 
oxia and low pH as inducing genomic instability by DNA repair activity being reduced 
[76], and enrichment for mismatch-repair deficient cells [77]. 

Thus, convincing evidence exists that a hypoxic and acidic environment will lead to 
genomic instability due to impairment of DNA repair processes. 

In addition, Suresh [78] has reviewed the biomechanics of normal and cancer cells 
and shows that transformed cells routinely have altered deformability. This effect is 
likely due to the differences in the actin microfilaments - the major structural element 
of a cell [79,80]. 

The cytoskeleton also consists of microtubules that have a bending stiffness 2.6 
X 10" 23 N m 2 which is about 1000 times stiffer than actin filaments [78]. Microtu- 
bules, unless stabilized by ligands or microtubule-associated protein, are in a con- 
stant dynamic instability process being polymerized and depolymerized with a half- 
life of about 2 minutes [81]. This dynamic instability phenomenon has been 
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mathematically modeled in detail by Sept et al [82], and Bolterauer et al [83] 
among others. The ends of actin filaments and microtubules are caped with ATP 
and guanosine triphosphate (GTP), respectively. GTP is created by the citric acid 
cycle in the mitochondria. Obviously, a mitochondria failure may reduce the con- 
centration of ATP and GTP, if other systems do not compensate for the mitochon- 
drial failure. A reduction in GTP results in a higher rate of depolymerization of 
microtubules [84] and a decrease in ATP concentration reduces the rate of growth 
of actin filaments [85]. Too few microtubules and/or microtubules being too short 
to participate in proper spindle pole formation can lead to mitotic catastrophe 
and/or potentially lead to chromosome instability. Bakhoum et al. [86] and Thompson 
et al. [87] describe the mechanism for this instability. Persistent mis-oriented attachment 
of chromosomes to the spindle microtubules leads to severe chromosome segregation 
defects. 

Conclusions 

Summary 

We have described an integrated system model for the progression of a healthy cell to 
a cancer state and some of the implications. The potentially aberrant state of the cell 
may start by an excess glucose or other nutrient external to the cell impacting the cell 
or by internal defects leading to metabolic enzyme redistribution processes. This excess 
nutrient is essentially a chemical potential difference between inside and outside of the 
cell creating stress. Through a process analogous to Rayleigh-Benard convection, stable 
molecular oscillators accumulate in the cytoplasm to exploit this chemical gradient. 
The continued activity of these oscillators results in mitochondrial destabilization, 
which may occur as a non-equilibrium phase transition. Once the mitochondria begin 
to perform aberrantly there will be a chemical imbalance in key components for micro- 
tubule assembly/disassembly. This imbalance is driven by Le Chateliers principle. The 
disruption in microtubule lengths and/or microtubule count will lead to chromosomal 
instability via kinetochore-microtubule dynamics finally leading to mitotic failure. In 
unlucky cases this will result in chromosome mis-segregation and cancer if mitotic ca- 
tastrophe does not occur. 

Clinical implications 

The above synthesis of ideas surrounding the subject of Warburg initiation, or the tran- 
sition from aerobic glycolysis to anaerobic suggests not only an avenue for treatment 
but also an avenue for prevention of cancer. 

We hypothesize that excess of glucose and sugar-like energy sources or metabolic en- 
zyme abnormalities, through a non-equilibrium phase transition (a symmetry breaking 
phenomenon) analogous to the Rayleigh-Benard convection, may cause a cell to prefer to 
process this energy source using substrate glycolysis. Continued excess substrate glycolysis 
will cause further phase transitions to disrupt the mitochondria via depolarization and also 
disrupt microtubule dynamics. When a cell then passes through mitosis, the chance of mi- 
totic failure is increased. When a cell enters mitotic failure, it may undergo an aneuploidy 
event [88,89]. All this suggests that a low glycemic diet would lower the incidence of can- 
cer, and may suggest a mechanism why metformin, which lowers blood glucose levels, is 
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associated with improved outcomes in diabetic cancer patients [90,91] and reduced risk of 
pancreatic cancer [92]. 

The above synthesis of ideas also supports targeting cells that have made the glyco- 
lytic switch. For instance, the work of Pedersen [93] and his colleagues (Ko et al. 
[69,94]; Mathupala et al. [70]) have used 3-bromopyruvate to inhibit glyceraldehyde 3- 
phosphate dehydrogenase (GAPDH), which effectively inhibits glycolysis [95]. In 
addition, 3-bromopyruvate may force, via Le Chatelier s principle, some reverse reac- 
tions to essentially deprive the cancer cell of substrate-created ATP. This leaves the cell 
little choice except to enter apoptosis. We further hypothesize that a Br derivative of 3- 
phosphoglycerate would similarly, though perhaps not as energetically, and perhaps not 
as toxically, facilitate via Le Chatelier s principle, a reverse reaction to deprive a cancer 
cell of ATP. 

Further, since microtubule dynamics are dysregulated by glucose oscillations and its 
associated pH oscillations, we speculate that metronomic dosing of microtubule poi- 
sons (e.g. nocodazol, taxol, vinblastine) would be an effective treatment strategy for 
cancers. 

In our list of "prescriptions" above many of these are already known or in use. These 
current practices are essentially "prediction" of our theory and were strictly based on 
biophysics with little detailed biochemical or cellular biochemistry being considered. 
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